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Abstract 

3d lattice studies have recently attracted a lot of attention, especially in connec- 
tion with finite temperature field theories. One ingredient in these studies is a 
perturbative computation of the 2-loop lattice counterterms, which are exact in 
the continuum limit. We extend previous such results to SU(iV) gauge theories 
with Higgs fields in the fundamental and adjoint representations. The fundamental 
SU(3)xSU(2) case might be relevant for the electroweak phase transition in the 
MSSM, and the adjoint case for the GUT phase transition and for QCD in the high 
temperature phase. We also revisit the standard SU(2)xU(l) and U(l) theories. 



1 Introduction 



Finite temperature gauge theories exhibit many interesting and important 
phenomena, such as the EW (electroweak) and QCD phase transitions. Un- 
fortunately these are not easy questions to answer, due in part to the infrared 
problem at finite temperature. However, recently a successful approach has 
been developed at least for the thermodynamical aspects of the EW phase 
transition (for a review, see [1]). This approach consists of a perturbative di- 
mensional reduction into a 3d effective theory [2-7] and of non-perturbative 
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lattice simulations in that theory [8-13]. Other cases where similar 3d sim- 
ulations are needed are QCD in the high temperature phase [14] and the 
Ginzburg-Landau model of superconductivity [15,16]. 

One ingredient in the non-perturbative 3d studies is the relation of the lattice 
and continuum regularization schemes [17,18]. This problem is analogous to 
the determination of Aj^/Aiatt in 4d [19]. In this paper, we extend previous 
3d results to theories relevant for extensions of the Standard Model. 

To be more specific, consider that one has to convert the results of lattice 
simulations to continuum observables. In 4d this is usually done by measuring 
some experimentally known observables and using them to fix the parameters. 
However, in 3d this is not possible since no such observables exist. Thus one 
has to use lattice perturbation theory [20] to explicitly renormalize the theory 
and find the relationship to the MS scheme observables in this way. For super- 
renormalizable theories the result will be exact in the continuum limit at a 
finite order in the loop expansion; in 3d one needs to go to 2 loops, and only 
the mass parameter gets renormalized. 

The task of relating lattice and continuum observables for 3d SU(2)-|-Higgs 
theories was previously discussed in [17,18]. In [18] the calculations were shown 
in some detail and results were given for SU(2) and U(l) theories with Higgs 
fields in adjoint and fundamental representations. Here we go to SU(A^) gauge 
fields with both types of scalar fields. 1-loop C(a)-improvement in 3d has been 
discussed in [21]. 

It should be noted that for some observables the 2-loop mass counterterms are 
not needed. Such case is, for instance, the determination of the discontinuity 
in (0V)/5'3 the phase transition point (m^) in a theory with only one 
mass parameter, if one does not need to know the numerical value of ml [21]. 
However, in most cases the value of is needed: one needs it in order to 
make a systematic comparison with 3d perturbation theory (the lattice and 
perturbative results should be compared for the same parameters), one needs 
it if there are several mass parameters in the theory as in the SU(3) xSU(2) 
model discussed in Sec. 5, and one needs it when connection to 4d physics is 
made, especially in theories sensitive to the value of T/A^^ such as QCD in 
the high-temperature phase [14]. 

This paper is organized as follows. In Sec. 2 we briefiy discuss some general 
aspects of the problem. In Sec. 3 we calculate the 2-loop lattice counterterms 
for the SU(A'")-|-adjoint Higgs theory and in Sec. 4 for the SU(A^)-|- fundamental 
Higgs theory. Some other fundamental theories, namely SU(3) xSU(2) with 
two Higgses and SU(2) xU(l), U(l) with one Higgs, are discussed in Sec. 5. In 
Appendix A we discuss the Feynman rules, in App. B some lattice integrals, 
in App. C the isospin contractions, in App. D the numerical values of the 
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constants appearing, and in App. E a unified representation for the g^,g^- 
parts of the counterterms. 



2 Calculations 



In order to find the relationship between the lattice and MS renormalization 
schemes we have to calculate some physical quantity in both schemes and 
equate them. As discussed in [18], the simplest choice is the minimum value 
of the effective potential. Practically this means requiring that the 0^-parts of 
the effective potentials coincide. 

In addition to the 0^-term, one needs to compute the mass-dependent part 
of the vacuum counterterm. Its definition is that if it is added to the lattice 
Lagrangian, >Ciat — > >Ciat + SV, then the mass-dependent part of the vacuum 
energy density agrees with the MS- result. The significance of the vacuum 
counterterm is that if £iatt = m^^V + ■ ■ then 



Hence, the relation of measured in continuum and on the lattice without 
adding SV to the Lagrangian is 

(0V)cont = (0V)latt + (2) 



This equation can be used to infer the non-perturbative value of (0V)cont from 
lattice simulations. Like for the mass counterterm, the 2-loop result is exact 
in the continuum limit. If one has a phase transition and is only interested in 
the discontinuity of then the vacuum counterterm does not contribute. 

When replacing SU(2) with SU(A^), the Feynman rules discussed in [18] get 

modified and obtain a more complicated form. In most cases this means only 
replacing e"^^*^ with f^^'", but sometimes more structure is added. The most 
crucial difference is that the symmetric tensors d'^^'" are no longer zero. The 
Feynman rules for the gauge and ghost-gauge vertices can be read from [20], 
with the replacements l/Aa'^ N/12a in Eq. (14.39), and {2/3){S^^S^^ + 
. . .) ^ {2/N){6^^6^^ + . . .) in Eq. (14.44). The sign of the ccAA-veitex in 
[20] should also be reversed: the ghost-gluon part of the Lagrangian is 
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where momentum conservation and due integrations are implied, and pi — 
I sin 7j = cos |rj. The other vertices are given in Appendix A. 

We will perform the calculations in the gauge with ^ = 1 in order to 
make the gauge propagator simple - the result is gauge fixing independent. 
As mentioned, we only need to go to 2-loop level to obtain the exact result. In 
[18] the results for most of the necessary integrals were given. In the present 
case, different components of the fields get different masses in the symmetry 
breakdown and we have to consider the corresponding integrals. However, it 
is easy to express all these new integrals using the ones defined in [18]; the 
most important cases are for completeness given in Appendix B. 



3 Adjoint Higgs 

Let us consider first a system with an SU(A'") gauge field and a Higgs scalar 
field in the adjoint representation of the gauge group. The Lagrangian on the 
lattice is 

« y i,j 

WTr$2 + Al(Tr$2)2 + A2Tr$^ (4) 
where Ui{x) — exp[iagAi{x)], Ai — AfT^, x + i = x + ae^, 

Pi,{x) = U,{x)U,{x + z)Ur\x + j)U-\x), (5) 

and a is the lattice constant. 

When a — > 0, this becomes the usual continuum Lagrangian 

Ccont = l^rFijFij + TrA$A$ + m^Tr^^ + Ai(Tr$2)2 + A2Tr$^ (6) 

where Fij = diAj - djAi + ig[Ai, Aj] and = + igiAi, 

We can write $ = Y1a=i^ $^T^, where the generators of the Lie algebra 
of the group SU(iV) are chosen to be Hermitian N x N matrices and are 



4 



normalized as TrT'^T^ — ^S"^^. If is a real A'"-dimensional vector, then 



T^'-'^ , ' (7) 
^2N{N - 1) V 



is Hermitian, traceless and properly normalized and can therefore be chosen 
to be one of the generators. 

We wish to relate the lattice observables to the MS scheme ones using the 
method given in [17,18]. In order to calculate the effective potential, we break 
the symmetry by shifting the Higgs field $ — > $ + vT^^~^. The effective 
potential in the adjoint case depends on the direction of the shift. However, we 
are only interested in the quadratic part of the potential which depends only on 
the magnitude v and thus any direction can be chosen. The Higgs-gauge cross 
term is cancelled by introducing an gauge fixing term — F^F^ /2^a^ ^ 
where 

F\x) = E [At{x) - At{x - i)] + ea^^r'^'^^-^$^(x). (8) 



If we define the projection operators 



pAB ^ ^AB _ pAB _ pAB. pAA ^ _ 

P^^ = -4 ^. 7^ . 7 + 2- — — P2 = 2(Ar - 1), 



we can write the propagators in the broken phase as 



pAB tdAB pAB 



k? -\- m\ k?' + ml k'^ + m^J^ 

/ pAB pAB pAB X /X 

where (...) is the same expression as in the gauge propagator and the masses 
are 



N 



2{N 



2 2 
:9 V , 



mi 



■m'^ + Ai + 



A^(A^- 1; 



-A2 
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iV2 - 3iV + 3 
Ar(Ar - 1) 



m^ = m^ + 3 Ai+ — ) u^. (11) 



It is now straightforward to calculate the renormalization countcrterms and 
to relate them to the ones obtained in the MS scheme. Let m^(/x) be the 
renormalized mass in the MS scheme and — m^{^) + Sm^di) the bare 
mass. The diagrams needed to calculate Sm^ are the same as in [18]. The 
Feynman rules for the vertices are given in Appendix A. The calculation of 
the isospin factors in the broken phase is discussed in Appendix C, and some 
typical lattice integrals are in Appendix B. For some diagrams, products of 
four structure constants are needed. Their values are given in Appendix A 
of [7]. 

Adding the contributions from the different diagrams together we get 



5m^ = - \2Ng^ + {N^ + l)Ai + {2N' 
1 



3) 



A, 



167r2 



N J Ana 
A. 



2Ng^ \^{N^ + l)Ai + (2N^ ~ ^) ]^ J ~ ^(A^^ + l)Ai 
N A^ 
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In— + C 

, afx , 



+2Ng'^ {N^ + l)Ai + {2N^ - 3) 



A2^ 

N 



T?-5 
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G"3^) ^S-4(5 + p) + 2/.i-/.4]} 



(12) 



where the constants C, 5, p, and E have been defined in [17,18]. Their 

numerical values are given in Appendix D. The n dependences of 5m^ and 
cancel as they should, which can be seen by comparison with Eq. (48) of [7]. 
As for SU(2), the finite parts of Sa/e arising from dimensional regularization 
in the MS case, cancel between the diagrams (vvv) and (vvs). The 2-loop 
1/a-terms cancel against contributions from the mass counterterms. 

The mass-dependent part of the vacuum counterterm is also straightforward 
to calculate, and it implies that 



Tr$2\ =/Tr$2\ _ (TV^ - 1)_ 

/ cont \ / latt ^ ' Stto 

-T^^(iV^-l)(>n^ + C+?-*). (13) 
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As a concrete example, consider SU(5). Using the numerical values of the 
constants E, S, p, (, Ki, K4 given in Appendix D, one gets from eq. (12), 



— w?'{ijl) 



47 

10(^2 + 26Ai + — A2 , ^ 

5 / 47ra 



^]_{[lO/(26A. + fA,)-52A;-15?A.A.-||A^ 



X (^In — - + c) + 124.01^^ + 5.7947^2 (^26Ai + y |. 
The combination appearing on the last row in eq. (13) is in eq. (D.5). 



(14) 



4 Fundamental Higgs 

For the fundamental theory, the lattice Lagrangian is 
+ \y\<t>\^)(t>{.x) - Re(t>\x)Ui{x)(t>{x + i) 

+ m^(/)V + A((/)V)^- (15) 
The corresponding continuum theory is 

Aont = ^TrF,,i^,,- + (A0)^(A0) + mVV + A(0V)', (16) 
where F^j = and A = 9^ + i^T^^f . 

In contrary to the adjoint case, the effective potential can be calculated in the 
fundamental case for an arbitrary shift (see, e.g., [22]). We make a shift such 
that 

= y- (17) 

For simplicity, we shall assume to be real (this assumption enters in the 
scalar propagators). The gauge fixing term is 

= E[^f (^) - - ^)] + - (18) 
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As mentioned, in practical calculations we use ^ = 1 to simplify the momentum 
algebra. The shifted field (p = d> — d> is written as 0„ = + iUa)/V^, 
a = 1, . . . , N. To get the propagators we define two more projection operators 
in addition to those in eq. (9), 



Ta/3 — Sal3 TTY^ ^'^^ ~ U "^"^ — ^ ~ — 1- (19) 

0T0 



Then the propagators are 



/ \ TaB , Lad I , TaQ LaB 
\UaUfj) — h , {uJa^^is) = h 3—, 

+ k"^ + m\ k'^ + k^ + 

/ pAB pAB pAB X /X 

Here c^, are the ghost fields, the expression (...) is the same as for the 
gauge fields, and the non- vanishing masses are 



„ 1 „ „ _2 2(A^ - 1) „ 
mf. —VTb' + 3Av^. 



f211 



The isospin algebra appearing in the graphs is for the most part the same as 
in continuum, and can be read from [22]. There are a few additional terms 
which do not appear in continuum since the vertex is proportional to powers of 
. Such terms come from the ccAA- vertex and the two AAAA- vertices. The 
isospin algebra of (ccAA) can be deduced from the graph (vv) in eq. (B.46) 
in [22] . The algebra related to the AAAA- vertex in eq. (A. 4) is 



is^B.CD ^ .Ac.BD ^ .AD.BC\ <t>^T^T^T^T^ 4> _ (A^^ - 1) (2Ar2 - 3) 
yo ±0 ±0 ) 

where we used that the masses in the propagators can be put to zero here. 
The algebra related to the symmetric a^-part of the four-gluon vertex coming 
from the plaquette is 



AT 



^AB^CD ^ ^abe^cde\ ^(^B^C) + iB^D) 



t2 



±(2N^ - Z),C''(ACA-U.^ ^ ^'^ -^>P^ -^> ,W (22) 
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Here we used eq. (B.9) according to which the result of the momentum integral 
is symmetric and quadratic in the masses in the continuum limit. 

After the momentum integrals, one can verify the cancellation of 1/a-terms 
against 1-loop counterterm contributions. The result left is 



N 



Attq 



X f In — + c) + 2Xg' 



4iV2 
,N+1 (\ 



N 



N-1 



N 



+9' 



47V2 



AN^ - 1. 



-E2 + 



-AN{N + l)p - 2{3N'^ - 1)5 + 2N\2ki - k^) 



(23) 



In fact, the 1-loop terms proportional to g'^ and the 2-loop terms proportional 
to g"^ in eqs. (12), (23) can be written in a unified way, and for completeness 
this formula is given in Appendix E. 



The relation of measured in continuum and on the lattice is 



N12 

(0V)cont = (0V)latt - -; 

47ra 



1 ,'(iV^-l)(.nA + c+iE=-4 (24) 



(47r)2 V a/j, 

Corrections of order a to the discontinuity of discussed in [21]. 



5 Other fundamental Higgs theories 



5.1 SU(N)xSU(2) 



It has been argued in [6,22] that an SU(3) xSU(2) gauge theory with two scalar 
fields might be relevant for the electroweak phase transition in MSSM. Let us 
here consider for generality SU(A'") xSU(2). The continuum Lagrangian is 



+ {DIU)\DIU) + mlU^U + Xu{U^Uf + -/H^HU^U. (25) 
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Here Df ^ di + igwf'A'^ and Df = di + igsT^C^ are the SU(2) and SU{N) 
covariant derivatives (t" = 0""/2 where a" are the Pauh matrices), gw and gs 
are the corresponding gauge couphngs, H is the Higgs doublet and U is the 
right-handed stop field. Since the two scalar fields are coupled only through 
the simple scalar vertex proportional to 7, discretization follows immediately 
from Sec. 4. 

The results for the counterterms can also mostly be read from Sec. 4. The 
reason is that when one makes a shift in H, say, then in the SU(iV)-sector 
one only needs to replace mfj — > + ^WH in the masses in eq. (21). 
Thus the result can be deduced from the mass-dependent part of the vacuum 
counterterm. The only new graph is the scalar figure-8 graph proportional 
to 7^. 

The results for the bare lattice mass parameters are 



ml = ml{^) - 



N 



-gl + 2{N + l)Xu + 2^ 



Ana 



167r2 



{N^ - 1) 



^ 4iV2 _iv + 3 

3 s 



, oiV+1 



Ar2 



^(37^^-4(Ar+l)A^-27^)] (in ^ + C 



N + 1 
N 



+ 



iV2 - 1 



/i/vv 



+4- 



1 


■4iV2 _ I 


47V2 


[ 4 



H ttE + + 1 



-AN{N + l)p - 2(3iV2 - 1)5 + 2Ar2(2Ki - k^) 



Attq 



51 
16 



gtv + 9A^5^ - 12A2, + {N^ - l)^gl - N^' 



6 



X In — + C + 
\ an ) 

3 



9Aif5^ + (iV'- 1)7^1 



+^^16 



15 4 

— E^ + -vrE + 5 - 24p - 22b + 8(2ki - K4) 



(26) 



(27) 



Using the numerical values in Appendix D, one gets for A'^ = 3, 



ml^mfjip) - 
1 



;9l 



Ana 



167r2 



+ 8\u + 27 
64 

+ ^At;g-16A^ + 375^-272 



InA + c^ 
ap J 
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+19.633^^ + 12.362A[/g + 1.73847^^|, 
mjj = mjjiii) - (^^gl, + 6Xh + 37^ 



(28) 



167r2 



47ra 

^9w + ^^H9^w - 12A^ + 87^ - 37^ 



a// 



+4.9941^^ + 5.2153A,j^^ + 4.635875||. 



(29) 



The vacuum counterterms do not get any corrections from 7, so that 



= (//<//),.., - A _ ^3A(l„ A + c + _ ^) . (30) 

The numerical value of the constant C, + S^/4 — 5 is in eq. (D.5). 



5.2 SU(2)xU(l) 



For the s^'^-part of standard electroweak SU(2)xU(l) theory, only the nu- 
merical values of the lattice counterterms have been previously given in lit- 
erature [9]. We give here the exact result in terms of the lattice constants 

The SU(2)xU(l) lattice Lagrangian is 



^iatt = -^ETr[l-P.,(:r) 



+|^E[i-pr(-)] 



-ia.i{x)(T-i 



+ m^^Tr$'f(a;)<l'(a;) + A iTr<l'^(x)<l'(a;) 

2 - 2 



(31) 



Here Pij is the U(l) plaquette with the U(l) link Ui{x) = exp[iai{x)], where 
ai{x) = ag'Bi{x)/2. Any positive number 7 in cq. (31) gives the same naive 
continuum limit. The Higgs field has been written as a matrix $ = (00), where 
is the standard Higgs doublet and = ia^(j)*. The continuum Lagrangian 
corresponding to eq. (31) is 



L = + ^BijB.j + {D,cf>yD,cl> + m'<f>^<f> + A(0t0)^ (32) 
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where Di — di + igAi + ig'Bi/2 and Bij — diBj — djBi. 

The calculation proceeds as before. The bare lattice mass parameter is 



Svra 



(3^^ + + 12A) 



1 



167r2 



2^ V32 8 6 



ttS 11 



S — 3p + 2k,i — 



+ 3A(3^^ + ^ 



2 , ^'2^/^__^ 



(33) 



The numerical value of the constant 2-loop part on the last three rows is 



167r2 



A 9 ,9 / 1.66290n m 

4.9941^^ - 0.88600/^'^ + (0.00932 + )g"^ 
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+5.2153A^2 + 1.7384A^ 



The relation of the order parameters is 



(34) 



(0V)cont = (0V)latt - 



27ra 



IGtt^ ^ V a// 

The SU(2)+Higgs theory follows as the limit g' 0. 



(35) 



5.3 U(l) 



For the U(l)+Higgs theory we extend previous results [18] by introducing 
an extra parameter 7 as in Ref. [9] and Sec. 5.2, allowing a simultaneous 
representation of the compact and non-compact lattice actions. 

The U(l)-|-Higgs theory is defined by 



2a4e2 



X) 



12 



" i 

where Ui{x) = exp[iaeAi{x)]. The continuum Lagrangian is 



(36) 



(37) 



with Di = di + ieAi. The bare mass parameter is 



m^^m^in) - (2e2 + 4A) 



47ra 



167r2 



(is^ - 1 - 25 - 4p + 1^ ^) + SAe^ f -E^ - 5 
where the numerical value of the last row is 



(38) 



167r2 



(_1.1068 + ^^)e^ + 4.6358Ae^ 
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The relation of the order parameters is 



(39) 



(40) 



The non-compact case follows as the limit 7 — > oo. 



6 Conclusions 



We have applied the previously developed method of deriving 2-loop lattice- 
continuum relations to 3d theories relevant for QCD and some extensions of 
the Standard Model at finite temperature. The main results arc the expres- 
sions for the bare lattice mass parameters in eqs. (12), (26), (27), with numerical 
values as given in eqs. (14), (28), (29). When an extrapolation to the contin- 
uum limit is made in the lattice simulations, these equations allow to extract 
the corresponding finite fixed value of the MS mass parameter m'^{fj,). Then 
one can compare lattice simulations with 3d perturbation theory, or if rri^^n) 
has been computed with methods of dimensional reduction from 4d, one can 
extract the physical temperature to which the simulations correspond. 
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A Feynman rules 



In this appendix we give the Feynman rules for the interactions of the Higgs 
field in the symmetric phase. The ones related to gauge and ghost-gauge ver- 
tices were discussed in Sec. 2. 

In the adjoint case, the vertices can be read from the action 



4 \ N J 8 

+ ^S{p + q + r + s)ip-^q)J^^'''f^^^^^ip)^^iq)Af{r)A^{s) 

_g2^jA,D,N^-ljBCD^^B^C (^.1) 

In the broken phase, the shift generates two new three-point vertices $$$ and 
^AA from those in eq. (A.l), but there is also an additional four-gluon vertex 



Sj, = -^^a^g^^f^^^f^'^'^'-^f^^f'^^^'^'-^AfAfA'^A^. (A.2) 



In the fundamental case, correspondingly, 

+ + g + r){pTq\(t>\p)T^(t>{q)Af{r) 

+ y + g + r + {p + q\ <t>^ {p)T^^T^mq)Af{r)Af{s) 
+ 5 Vc^ (^^'^T^T^ct) + 0tr^T^0) . (A.3) 
In the broken phase, the shift generates again the extra vertices 000, 074^4, 
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and there is the additional four-gluon vertex 



D 



(A.4) 



B Lattice integrals 



In [17,18], the lattice integrals were worked out for mass combinations specific 
to SU(2)+fundamental Higgs. Let us here describe how these generalize to the 
present case. Let 



Sij 



FijkiP-, q, r) = Sik qk{pj - Tj) + SkjPjivi - Qi) + 6jiri{qk - Pk), 



(B.l) 
(B.2) 



where the superscript in V^] is to remind that we are in the Feynman gauge. 
Then there are contributions from the following integrals [22] : 



D'yyy = j dp dqVli{p\ Ml)V}^{q\ Ml)VL{r\ Ml) 

X Fijkip. q, r)Fimn{p, q, r) 



1 



(47r): 



■E 

i=l 



-6M^ ( In — + + -Afi(87r - 15E) 
^ an J a 



TT, 



- + -E^ s + 65 + 6p - 4ki + 2^4 

2 4 3 



^vvs = jdpdqVli{p\Ml)V^^{q\Ml)S{r\m')A5,,n5irari 



12(ln— + C] + 2-E2 



L'^sv = / dpdqS{p\ ml)S{q\ m^)y.J(r^ M^){p - q),{p - q), 



1 



(M^ - 2m? - 2ml) (in — + C ) (^i + ^2) 

\ afj, J a 

3. 



(47r)2 

+ ^M(^7r - ^Sj + 25(m? + ml) + ApM^ 
D'ggv = jdpdqS{p\ml)S{q\ml)Vl^{T\M^){-p,pjq^(U) 

^Jl^\\M''-ml-ml)(\n — + C\+-MU--T. 



(B.3) 



(B.4) 



(B.5) 



(B.6) 
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where r = —p — q, m = m{l + ^am), ^ ~ 0.152859 (however, ^ cancels in the 
final result, see [17,18]), and 



n/a 



— 7r/o 



d^p 
(2^ 



(B.7) 



The contributions from the figure-8 graphs follow directly from the SU(2) 
case in [18], apart from the graph (vv). Let us note that there are two kinds of 
contributions from this graph, one proportional to f^'^^f^^^{A'^A^){A'^A^) 
(this is the isospin factor appearing in continuum) and the other to the factor 
in eq. (22). In the former case the integral appearing is 



Di(^)= jdpdqS{p\Ml)S{q\Ml 



(47r)2 



E 

i=l 



1 



28 



10 



-Mi ( -12E + —n] + Mfi-T.' - — ttS 



and in the latter case 



(B.8) 



D]^^ jdpdqS{p\Ml)S{q\Mi 



4~2~2 

--apq 



(B.9) 



C Isospin contractions 



When calculating the adjoint case mass countertcrm in eq. (12) one has to 
evaluate 2-loop Fcynman diagrams with complicated vertices. There are basi- 
cally two ways to evaluate the isospin factors related to the diagrams. Either 
one can use the Fiertz identity 

T^T^i = \ - ^) (C.l) 

and the projection operators in eq. (9), or one can write the fields in component 
form and calculate products of the structure constants of the Lie algebra. We 
give here the necessary products for the latter strategy, which usually leads 
to simpler algebra. 

Let [A < N"^ — 1) be Hermitian traceless N x N matrices with TxT'^T^ ~ 
They form a representation of the Lie algebra su(iV). Let us define 
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~^ according to eq. (7). The matrices can be chosen such that the 
projections Pf"^ separate the index set / = {1, . . . , N"^ — 1} to three subsets 
Ji, I2 and 73. The matrices corresponding to set Ji form a representation of 
su(A^ — 1) and Is contains only the index N"^ — 1. 

Since the Feynman rules for the vertices contain factors of the form f^^'" , 
jABE jCDE ^ABE^CDE obvious that One needs to calculate sums like 



A£U Beij Ceik 

^S=E E E^^^^^^^^> (c-3) 

A€li B€lj C€lk 

^J? = E E E rf-^^^rf^^^- (C.4) 

Aeh Beii Cei 



Using the known results 



jACDjBCD^^AB^^ ^AC D ^BC D ^ ^AB ^ ^AAB ^ ^^^ ^^ 

we can deduce the following values for Fijk, D^jl and D^j^: 



Fni^N{N-l){N-2), F,22 = N{N-2), F223 ^ N, (C.6) 
(1) ^ (iV + l)iV(7V-2)(7V-3) ^ N-2 
111 7V-1 ' iV-1' 

dSs^^-^, ^S = 2||^, (C.7) 



dS^2N^-^^, Df^^-2{N-2)\ 

2 



The results are symmetric in the permutation of the indices. All other sums 
are zero. 



D Numerical values 



Let us use the shorthand s = sin. Then the numerical constants appearing in 
the 2-loop calculation are as follows: 
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7r/2 

^ / £x — ^ = 3.175911535625, 



-tt/2 
7r/2 



^=^. d'^d^ ^^ 2 ^'/''^^ V 2 - 1-942130(1), 

— 7r/2 

''^/^ ''2 2 4 

Ei "5 (Xj -|- yij Ei -^j I 1 



— 7r/2 

= -0.313964(1), 

— 7r/2 
— 7r/2 

The constant E can be expressed in terms of the complete eUiptic integral of 
the first kind [17]. The number in parentheses in the other constants estimates 
the uncertainty in the last digit. 

In [18] two more constants K2 and K3 were defined. However, they appear in 
the calculations only as a sum K2 + Ks and can thus be eliminated using the 
relation 

2{5 + p) = K2 + K3 + K4. (D.2) 

This follows from the trigonometric identity 

sin^ X + sin*^ y + sin^(a; + y) + A sin^ x sin^ y sin^(a; + y) 
—2 ^sin^ xsin^ y + sin^ a; sin^ (a; + y) + sin^ y sin^(a; + y)^ = 0. (D-3) 

In addition, there is the constant ( for which one to our knowledge cannot 
write a direct expression in momentum space: 



( = lim 



1 ''f d^xd^y , 3 1 

' ^ - In - - 



47r4 J (Ei s^Xi + z^) Ei s^{xi + yi) Ei s^yi z 2 

— 7r/2 

= 0.08849(1). (D.4) 



For the numerical evaluation this integral can be reduced to a four- dimensional 
one [17]. Note also that for the combination appearing in the relation of lattice 
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and continuum order parameters in eqs. (13), (24), (30), (35), (40), one gets 

C + |E2- (5 = 0.66796(1). (D.5) 



Finally, let us point out that the results above follow from a straightforward 
numerical integration. If needed, more precise values could be obtained using 
the position space techniques developed in [23]. We these techniques one can, 
for instance, write a closed expression for (. Noting from eq. (11) in [24] that 
replacing the mass am — 2z in eq. (D.4) with an external momentum ap 
causes the numerical factor 1/2 to be replaced with 3/2, one gets 

C = - In 6 + E [{infGixf + ^Hix)] . (D.6) 



Here (we have put a = 1 as in [23]) 



The functions G{x), H{x) are related by H{x) = 2J2i G{x)-G{x-i) /(Ei Xi) 
for 7^ 0, and have at x — the values 

G{0) = ^, H{0) = 1.67338930297(1). (D.8) 
47r 



The sum in eq. (D.6) is convergent since the leading asymptotic behaviours 

G(x) H(x) ~ ^ (D.9) 



cancel in eq. (D.6). To evaluate the sum one needs an efficient method of 
calculating G{x) very precisely [23]. 



E A unified representation for the gauge peirt of the counterterms 



Let F^Q — —if^^^ be the generators of the adjoint representation, and G^ 
those of a general representation. We define T such that for a real (adjoint) 
representation, 

TrF^F^ = r5^^, (E.l) 
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and for a complex (fundamental) representation, 



(E.2) 



The quadratic Casimir C2 is defined by 

{G^G^)ab — C25ab- 



(E.3) 



In addition, the number of real scalar field components is denoted by d. Then 
for the fundamental representation, T = 1, C2 = (A^^ — 1)/2N and d — 2N, 
and for the adjoint representation, T — C2 — N, d — N"^ — 1. 

In the limit that the scalar self coupling constants vanish, one can then write 
the bare mass parameter in a unified way: 



(5m^ = -2/C2 



Ana 



9^ 



167r2 



Co 



+ \C2) + (y - 4 J + ^ - ^2 



-2{N + T)p-2{N + C2)5 + N {2ki - k^) 



The mass-dependent part of the vacuum counterterm is 



5V = -m^d 



+ TTT^Ciiln — + C + -r - ^ 



Stto {Any 



(E.4) 



(E.5) 
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